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Abstract 



^ The combination of density-functional theory with other approaches to the many-electron prob- 

Q lem through the separation of the electron-electron interaction into a short-range and a long-range 
contribution (range separation) is a successful strategy, which is raising more and more interest in 

g 

I recent years. We focus here on a range-separated method in which only the short-range correlation 

I ... ... . . . ... . ... ..... O.... ....... 

I We consider the paradigmatic case of the H2 molecule along the dissociation curve, finding encour- 

1—1 aging results. By means of very accurate variational wavefunctions, we also study how the effective 
> 

O electron-electron interaction appearing in the Overhauser model should be in order to yield the 

^ exact correlation energy for standard Kohn-Sham density functional theory. 
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I. INTRODUCTION 



Kohn-Sham (KS) density functional theory (DFT) (see, e.g., [T]) is a successful method 
for electronic structure calculations, thanks to its unique combination of low computational 
cost and reasonable accuracy. In the Kohn-Sham formalism, the total energy of a many- 
electron system in the external potential Vne = J2i "^nelrj) is rewritten as a functional of the 
one-electron density p(r), 

E[p] = Ts[p] + j dvvneir) p(r) + U[p] + (1) 

In Eq. ([T|, Ts[p] is the kinetic energy of a non-interacting system of fermions (usually called 
KS system) having the same one-electron density p of the physical, interacting, system. 
The Hartree energy U[p] is the classical repulsion energy, U[p] = ^ f dr f dr' p{r) p{r')\r — 
r'l^^, and the exchange-correlation functional Exc[p] must be approximated. Minimization 
of Eq. ([T]) with respect to the spin-orbitals forming the KS determinant lead to the KS 
equations. Thus, instead of the physical problem, in KS DFT we solve the hamiltonian of a 
model system of non-interacting fermions, and we recover the energy of the physical system 
via an approximate functional. 

Despite its success in scientific areas ranging from material science to biology, approx- 
imate KS DFT is far from being perfect, and many fundamental issues still need to be 
addressed. In particular, KS DFT encounters difficulties in handling near-degeneracy cor- 
relation effects (rearrangement of electrons within partially filled shells), and in taking into 
account long-range van der Waals interaction energies (crucial, e.g., for layered materials and 
bio molecules) . In principle, all the shortcomings of KS DFT come from our lack of knowl- 
edge of the exchange-correlation functional, and a huge effort is put nowadays in trying to 
improve the approximations for E^dp] (for recent reviews see, e.g., [21 E]). 

An alternative strategy to overcome the problems of DFT is range separation: the 
electron-electron interaction is split into a long-range and a short range part, and the two 
are treated at different levels of approximation [lElElElElElIinillllinillSlIIlIISlIini 
[I71[ISl[I2ll2niEIll221l23ll2ai2Sll2Sll2Zl. Prof. Hirao has been a pioneer in this field, 
investigating the effect of range separation on the exchange energy with remarkable success 
(see, e.g., HI El [18]). 

The variant of range separation that we consider here [6l [71 El El HDl El [121 HSl [HI [151 

[T6l [T7] can be viewed as a way to remove the constraint that the model system be non- 
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interacting: instead of the KS system, one can define a long-range-only-interacting system 
(whose wavefunction is thus multideterminantal) having the same density of the physical 
system. The remaining part of the energy is then approximated with a short-range exchange- 
correlation functional. The resulting long-range-only hamiltonian, being weakly interacting 
(and without the electron-electron cusp), can be treated at a reasonable computational cost 
with standard wavefunction methods: in general, the needed configuration space to achieve 
good accuracy is small, and often second-order perturbation theory suffices. At the same 
time, this long-range interaction, albeit small, can make the corresponding wavefunction 
capture near-degeneracy effects and long-range van der Waals energies. Provided that the 
energy functionals are correctly redefined, there is no double counting of the energy, and the 
method is in principle exact, as it is KS DPT. 

As mentioned, this range-separated multideterminant DPT needs an approximation for 
the short-range exchange-correlation functional. One can follow the same path as for KS 
DPT: start with the local-spin-density approximation (LSDA), consistently constructed as 
the difference between the standard LSD functional and the exchange-correlation energy 
of an electron gas with long-range-only interaction [28], and then add gradient corrections 
(GGA) [111 [121 El EH], and eventually meta-gradient corrections (mGGA). However, this 
path, which proved highly successful for KS DPT, may not be the best for a scheme in which 
long-range correlations are explicitly taken into account by wavefunction methods. Indeed, 
in most cases there is no improvement when passing from LSDA to GGA [HI |T2l EO], with 
the exception of hydrogen-bonded complexes [^ . 

In recent years we have extended the "Overhauser model", an approximate method to 
calculate the short-range part of the pair density in the uniform electron gas, to systems 
of nonuniform density [321 |33l [SH |35] , finding that it yields an accurate description of the 
short-range part of the spherically- and system-averaged pair density (intracule density) of 
small atoms. In Ref. [55] we have combined the Overhauser equations with the Kohn-Sham 
equations in a self-consistent way, recovering full CI total energies within 1 mH for the He 
isoelectronic series. In Ref. [3l] we have shown that, unlike all the available correlation 
functionals [36j, the model works equally well for the high-density limit of the He and 
the Hooke's atom series. Thus, on one hand, the Overhauser model seems to be a very 
good candidate to construct short-range correlation energy functionals. On the other hand, 
we have tested it only on systems dominated by dynamical correlation: in the He atom. 
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the Overhauser model yields essentially the exact KS correlation energy. However, when 
we move to systems with strong static correlation we expect the Overhauser model to be 
unable to yield good results. The combination of the Overhauser model with range-separated 
multideterminant DFT seems then natural: it can be viewed as a way to produce an adapted 
short-range correlation functional for the range-separated multideterminant DFT, or as a 
way to add the description of static correlation to the Overhauser model. 

In this work we combine the Overhauser model with range-separated multideterminant 
DFT, applying it to the case of the H2 molecule along the dissociation curve, thus analyzing 
also the case of strong static correlation, as the dissociation limit is approached. The paper 
is organized as follows. After briefly reviewing in Sees. [TT] and III the basic equations of 



range-separated multideterminant DFT and of the extended Overhauser model, we first 



analyze in Sec. IV, using very accurate variational wavefunctions [371 l38| [39] . how the 
"exact" electron-electron interaction which appears in the Overhauser model (and that it is 
usually approximated with a physically-motivated interaction) should be as the H2 molecule 
is stretched. This analysis shows the difficulty of modeling static correlation within the 
Overhauser model. Since the model is only able to describe correlation, we combine it 
with a generalized OEP scheme for multideterminant DFT, which is described in Sec. IVl 



The combination of the two methods is then presented in Sec. VI, with results for the H2 



molecule. The last Sec. VII is devoted to conclusions and perspectives. 



II. MULTIDETERMINANT DFT VIA RANGE SEPARATION 

Hohenberg and Kohn [ID] introduced a universal functional of the density F[p\, which 
can be written clS cl constrained minimum search 



F[p] =min(^|T + Ke|^). (2) 

In Eq. (2| the expectation of the kinetic energy operator T = Si Vf plus the Coulomb 
electron-electron repulsion operator V^e = X]i>j is minimized over all wavefunctions 

yielding the density p. The universality of the functional F[p\ stems from the fact that T 
and Vee are the same for every electronic system of given particle number N = J p{r)dr. 
Kohn and Sham [12] introduced another functional, Ts[p] of Eq. (Ilj), by replacing Vee in 
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T,[p] =min($|f|$), (3) 



Eq. g with zero, 



and used Ts[p\ for approximating an important part of F[p\. In Eq. (|3]), and in the rest of 
this paper, $ denotes a non-interacting wavefunction (thus in the majority of cases a single 
Slater determinant). Similarly, we can introduce a functional -Flr[p] for a long-range-only 
interaction W^j^ (here chosen using the error function, with the real parameter p governing 
the cutoff of the short-range part), 



by defining 



In this way we have 



Fl^[p]= m:m{^^\f + Wl^\^n. (5) 



hm Fl^[p] = F[p] (6) 

>oo 

lini^L^Rb] = T^.H- (7) 
We can then write the total energy of a given many-electron system as 

E[p] = Fi:^[p] + J drvUr) p{r) + j dv j dv' p{v)p{vf^^^^^^^ + i?,^Jp], (8) 

where erfc(x) = 1 — erf(x) is the complementary error function. As in KS DFT then, 
minimization is performed over the wavefunction \E'^, 

+ JdrJ drVM(r)p^.(rO '''^'^^^^^ ~ ""^'^ + Ei:,[p^.]j, (9) 

where p^M is the density corresponding to \E''^. Eq. (joj) yields an effective, long-range-only- 
interacting hamiltonian to be solved with a chosen wavefunction method. The short-range 
exchange-correlation functional E^^[p] is then defined as the energy needed to make Eq. ([s]) 

EM = F[P] - F{:k[p] - jdvj rfrV(r)p(rO ''^'^^^^[';^'^'^ (10) 
For instance, the correct LSD approximation to E^J\^p] is 

^x^c''°[p] = / P(r){e.e(PT(r),Pi(r)) -6^,(p,(r),p,(r))}, (11) 



where exc{P'\{r), pi{r)) is the exchange-correlation energy per electron of the standard 
uniform electron gas (with Coulomb electron-electron interaction) and e^^(p|(r), p|(r)) is 
the exchange- correlation energy per electron of a uniform electron gas with interaction 
erf(/iri2)/ri2 [28]. 

An exact expression for E^^[p] is found from the adiabatic connection formula fl6[ E3] : 

E^M = f dp' 4rrrlJ^\ru)^e-^''^'i^dn,- I dr J rfrV(r)p(rO^^^^^^^^, 

(12) 

where f^{r 12) is the spherically and system- averaged pair density (intracule density) obtained 
by integrating \^}^^\'^ over all variables but ri2 = |r2 — ri|, 

r(ri2) = y"|^''(ri2,R,r3,...,r^)r^rfMr3...rfr^, (13) 

(T1...(7JV 



with R = (ri +r2)/2. The gaussian damping appearing in Eq. ( 12 ) comes from the derivative 
of the long-range interaction erf(/iri2)/ri2 with respect to p, and shows that the exchange- 
correlation energy is determined by the short-range part of the intracule density. Notice that 



when yU = Eqs. (12) yields the KS exchange-correlation energy functional from a nonlinear 



adiabatic connection IT6l H3i. 



III. THE EXTENDED OVERHAUSER MODEL 

The extended Overhuaser model consists in writing an effective Schrodinger-like equation 
for the intracule density f{r 12) of a given system. The basic idea is the following [32l 
ESI ES]- We start with the observation that the intracule density f{r 12) couples to any 
electron-electron interaction operator depending only on the interelectronic distance, W = 
J2i>j '^(kj — '^jl)) the same way as the density p(r) couples to any local one-body potential 
operator V = i-^-, 

{^\W\^) = J dri2f{ri2)w{n2), (14) 
= [ drp{r)v{r). (15) 



We can then follow the Hohenberg and Kohn philosophy but with the roles of p{r) and 
f{^i2)i and of Vne and V^ej interchanged. That is, in analogy with Eq. we can define a 
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system-dependent functional G[f], 

Gi/] = mm(^|f + Ke|^), (16) 
so that the total energy of a given physical system is equal to 

E[f] = G[f] + [ dvj-^. (17) 

Like Kohn and Sham, we can define another functional by setting Vne equal to zero in 



Eq. (16) 



Tf[/]=min(vl/|f|M/). (18) 

The functional Tf[/] corresponds to the internal kinetic energy of a free (zero external po- 
tential) cluster of fermions having the same intracule density of the physical system. The 
fermions of this cluster interact with an effective interaction Wef[{r 12) which has the same 
role of the KS potential for the KS system. In practice, this effective interaction must be 
approximated. Moreover, for N > 2 electrons the cluster equation become a complicated 
many-body problem, so that other approximations are needed. As in the original Overhauser 
model for the uniform electron gas H5] , we can approximate the cluster equation with a 
set of radial geminals gi{r 12), 



1 i{i + 1] 



r o-ri2 H 2 ^ Wesiru) 

ri2 drf2 



gi{ri2) = eigi{ri2) 



J2miri2)\' = firu), (19) 

i 

whose occupancies i)i must be defined (e.g., in a determinantal-like way as in the original 
Overhauser model |15]). In practice, trying to solve the whole many-electron Schrodinger 



equation by means of Eqs. (17)-(19) is a daunting task. The idea is rather 



to couple this "average-pair-density-functional theory" with a density functional scheme: 



Eqs. (19) can be generalized to any f^ir^) along the adiabatic connection of DPT. In 



Refs. [321 [35J we started from the effective interaction w^^{ri2) which, when inserted in 



Eqs. (19), gives the intracule density corresponding to the Kohn-Sham system, /ks(?"i2) 
(that can be obtained from the KS determinant). We then wrote an approximation for 
Wgg(ri2) along the long-range adiabatic connection of DFT as 

<niru) = wfi{n,) + w:^^{n,). (20) 
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The only term that needs to be approximated is then wl^lr 12), an effective interaction that 
should essentially "tell" to the intracule density that, while the electron-electron interaction 
is turned on (i.e. as fi increases), the one-electron density p(r) does not change. As the 
information on p(r) has been "washed away" in the integration over the center of mass R 



of Eq. (13), this constraint can be imposed only in an approximate way. For two-electron 
atoms, for which Eq. (19) is exact with one geminal j33j, 9 = \/f, a simple approximation 
for <ff^(ri2) is [321 132 [35] 



w. 



eS ' 



ri2 



erf(/iri2) 

ri2 



3 ' 



-1 



erf(/i|ri2 - x|) 



fix. 



|x|<rs 



ri2 



(21) 



where fg is a screening length associated to the radius of a sphere containing on average one 
electron [HI 133 HE]. The physical idea behind Eq. (21) is to mimic the constraint of fixed 
one-electron density by screening the electron-electron interaction over a length associated 
to the "space" available to each electron (which is determined by the density). Indeed, 



for the He isoelectronic series Eqs. (12), (19) and (21), combined self-consistently with the 
Kohn-Sham equations, recover the full CI total energy within 1 mH [321 |3ll [35] . 



IV. THE OVERHAUSER MODEL FOR THE H2 MOLECULE: HOW THINGS 
SHOULD BE 



For a closed-shell physical electronic system (atom, molecule) with N = 2 particles, the 
Schrodinger equation describing the internal degrees of freedom of a cluster of fermions 
having the same intracule density f{r 12) is exactly given by [331 EH [35] 

1 (P 



ri2 dr"^ 



-ri2 + Wesiru) 



12 



ri2 



(22) 



As a first study, we calculate and analyze the "exact" Overhauser interaction Wes{f 12) at 
full coupling strength (i.e., for electron-electron interaction l/ri2, corresponding to /i = 00) 
for the H2 molecule at different values of the internuclear distance R, and we compare it 



with the approximation of Eq. (21). To this purpose, we need extremely accurate intracule 



densities /(ri2), which are described in the next Subsec. IV A 



8 




0.016 
0.012 
0.008 
0.004 




'12 



physical x 5 
KS 



R = 4.5 




10 



'12 



0.016 
0.012 
0.008 
0.004 




0.016 
0.012 
0.008 
0.004 




physical 
KS 




R = 3.0 



012345678 



'12 



physical x 10 
KS 




R = 6.0 



8 10 12 



'12 



FIG. 1: Intracule densities f{ri2) for the H2 molecule at different internuclear distances R for the 



physical system (from the accurate variational wavefunction described in Subsec.[lV A) and for the 



KS system (from the density corresponding to the same accurate variational wavefunctions) . 



Intracule densities from accurate variational wavefunctions 



We use the accurate variational wavefunctions of Refs. ESI ES], which are expanded 
in explicitly correlated gaussian geminals, 

K 

^(ri,r2) = (1 + A2)(l + y J]cfe^fe(ri,r2) (23) 



V'fc(ri,r2) 



k=l 

g-Ofc|ri-rAfcPg-A|r2-rflfe|2g- 



(24) 



where r^k and vsk are centers that lie on the internuclear axis, P12 nieans permutation 
of ri and r2, and ie is the inversion operator with respect to the center of the molecule. 



The parameters appearing in Eqs. (23)- (24) are determined variationally by minimizing the 



energy with the conjugate gradient method (for more details on the wavefunction and the 
algorithms employed, see Refs. 



[39]). The expansion length K = 1200 in Eq. (23) is 
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FIG. 2: The same intracule densities of Fig. [l] multiplied by the volume element 47rr^2- 

used, resulting in energies with the extraordinary accuracy of 10"^*^ Hartree. 

The intracule densities f{ri2) from these extremely accurate wavefunctions can be easily 
calculated, since all the needed integrals are analytic. We also calculated the one electron 
densities p(r), and the intracule densities /ks(^i2) corresponding to the KS system, which 



can be obtained by inserting in Eq. (13) the KS wavefunction ^y^p{ri)y^p(r2). In Fig. [ijwe 
show the intracule densities /(ri2) and /ks('"i2) for the internuclear distances R = 1.4, 3.0, 



4.5 and 6.0 a.u. Although mathematically the wave function of Eqs. (23)-(24) is cuspless. 



we see that the very elaborate ansatz permits to describe the exact linear behaviour of 
the intracule density for ri2 0, up to extremely short distances. Fig. [2] shows the same 
quantities multiplied by the volume element 47rr^2- This figure better visualizes the transition 
from dynamical to static correlation. In Fig. [3] we also report the same quantities in the 
extreme stretched case, R = 20, obtained from the simple Heitler-London wavefunction. 
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FIG. 3: The intracule density /(ri2) multiplied by the volume element 47rr^2 t^ie H2 molecule 
in the extreme stretched case R = 20. The physical /(ri2) has been calculated from the simple 
Heitler-London wavefunction, and /ks(^12) from its corresponding density. 
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FIG. 4: The "exact" effective Overhauser interaction U'efr(ri2) for the intracule densities of the 
physical and of the KS systems of Fig. [T] 



11 




0.5 1 



•12 




O 0) 




•12 




•12 



•12 



FIG. 5: The difference w^g{r 12) between tlie "exact" effective Overliauser interaction for the in- 
tracule density of the physical and of the KS systems of Fig. |4| The Coulomb repulsion 1 /ri2 is 
also reported. 
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FIG. 6: Same as Figs. |4]|5] for the extreme stretched molecule, using the simple Heitler-London 
wavefunction. 
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B. Accurate Overhauser potentials 



From the accurate intracule densities of the previous subsection we can calculate, by 
inversion, the corresponding "exact" Overhauser interaction Wesiru), 

Weff(ri2) = — (ri2VJ(ru)) + const. (25) 
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The inversion of Eq. (25) is done numerically, by finite differences. In Fig. 4 we report 



the effective Overhauser interactions that, when inserted in Eq. (22), give the physical and 
the KS intracule, corresponding, respectively, to /i = oo and /i = along the long-range 
adiabatic connection of Sec. [ll] (or to A = 1 and A = along the usual linear adiabatic 
connection in which Vee is simply multiplied by A). We see that Weff(ri2) for large ryi goes 
to the same constant for both the KS and the physical system, as it should be [35] (of course 
if we go to too large ryi we start to observe the wrong harmonic wall due to the gaussian 
asymptotic decay of our wavefunction). The difference between the effective Overhauser 



interaction for the physical and the KS system gives = ty^g of Eq. (20), and is 



reported in Fig. |5| where also the Coulomb repulsion l/ri2 is shown. From this figure, we 
see that, when the system is still dominated by dynamical correlation, as in the case R = 1.4 
and R = 3, Wgfj(ri2) is essentially a screened Coulomb interaction. That is, for short-range 
it behaves as 1/ti2, and then for large ri2 goes to zero much faster than l/ri2. In such cases. 



the approximation of Eq. (21), which at yU = oo reads 



<s{ri2) =4 + ^-4 ri2 < r, 

<ff(ri2)= ri2>r,. (26) 

can work reasonably well, with a screening length fg ~ R. However, as R grows and the 



system starts to be dominated by static correlation, we see that the approximation of Eq. (26 ) 
cannot work: the "exact" Wgfj(ri2) still decays much faster than l/ri2 for large ri2, but at 
short range is more repulsive than the Coulomb interaction! I.e., we need an "overscreened" 
interaction. This is completely evident in the extreme stretched case i? = 20 of Fig. [6} again 
obtained from the simple Heitler-London wavefunction. 
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V. GENERALIZED OPTIMIZED EFFECTIVE POTENTIAL METHOD FOR 
MULTIDETERMINANT DFT 



In recent years, the focus of a large part of the scientific community working on improving 
the approximations for Exc[p\ has shifted from seeking exphcit functionals of the density (hke 
the generahzed gradient approximations), to imphcit functionals, typically using the exact 
exchange Ey\p\, which is only explicitly known in terms of the Kohn-Sham orbitals 0i(r). The 
corresponding Kohn-Sham potential must then be computed with the optimized effective 
potential (OEP) method (for a recent review, see [H]). The OEP scheme can be generalized 
to the multideterminant range-separated DFT by first noticing that we can divide E!^J\^p\ 
into exchange and correlation in two different ways [10]: we can define the exchange energy 
with respect to the KS determinant $, 

= mVee - Wl^\^) - jdv j ^rV(r)p(rO ''^'^^^^'~^'^'\ (27) 

and then define the usual correlation energy functional -E^fp] as the energy missed by the 
KS wavef unction, 

mp] = E'.Ap\-mpi (28) 

but we can also define a multideterminantal (md) exchange functional [TO] by using the 
wavefunction 

KUp\ = (^l^ee - K^l^n -jdvj drV(r)p(r')^^^^^^^^, (29) 

and then a corresponding correlation energy that recovers the energy missed by \E'^ (which 
is smaller than the energy missed by the KS determinant $), 

Kr.M=E'lAp\-ElM- (30) 

Then, with this latter definition of the correlation energy, the generalized OEP-like scheme 
for multideterminant DFT becomes [TU] 

= inf {(v|>(;.|f + te + Kel^C.) + i^c%[P<j} , (31) 

where is obtained by solving the Schrodinger equation corresponding to the hamiltonian 

H^" = f + W^^ + V^", = v^'iv,). (32) 

i 
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Notice that this multideterminant OEP scheme is different from the one recently proposed 
In Eq. (31) the weak long-range interaction Wj^j^ automatically selects the 



in Ref. 



configuration space needed to yield an accurate solution for the hamiltonian H^^ of Eq. (32), 
while in Ref. [H] the configuration space is chosen essentially by hand, using physical and 
chemical intuition. 



In the next Sec. VI we use the Overhauser model to approximate E^^^^ip], and we apply 



our combined formalism to the case of the H2 molecule. 

VI. MULTIDETERMINANT DFT COMBINED WITH THE OVERHAUSER 
MODEL 



From the adiabatic connection formalism we can easily write an exact formula for -E^ mdH 

2 



djj! 



Airr 



12 



r {ri2)-nru) 



e ^'^ '^^ dr 



TC 



12, 



(33) 



which shows that El!^^[p] is determined by the change in the short-range part of the in- 
tracule density when the electron-electron interaction increases from erf (/iri2)/ri2 to the full 
Coulomb repulsion l/ri2. By adding and subtracting /ks(^i2)? Eq- ( [33| ) can also be written 

as 



KnJp] = mp] 



4^^?2 [r(n2) - /Ks(ri2)] 



erfc(yur 



12 



dr 



ri2 



12, 



(34) 



where Ej^[p] is the correlation energy of Eq. (28), defined with respect to the KS determinant. 



We computed f^{ri2) within the Overhauser model, Eq. (19) with one geminal g = \fj . 



using the simple screened potential of Eqs. (20)-(21) with the screening length = R. For 



each internuclear distance i?, the intracules f^ir^) have been calculated for 33 values of p 
between p = 0.01 and p = 20. By numerical integration we then computed 



dE^[p] 
dp 



^Trrl, [r(ri2) - /Ks(ri2)] ^e"^ dr^, 

V71" 



(35) 



and we fitted the values of ^^^M ^ith the derivative of the function 

Ofl 



mp] 



04 



(36) 



(1 + 62^,2)5 

which has the correct asymptotic behaviors [^. We also computed, again by numerical 



integration, the second term on the right-hand-side of Eq. (34) in order to obtain -^^^^^[p]- 
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We then implemented the generahzed OEP scheme of Eq. (31) by first minimizing the 
effective potential t>^(r) at the "generalized-exchange" -only level, and by adding E^^^[p] 
only as a final correction. Since E^y^dip] ^^^Y small, we do not expect substantial changes 
by implementing a full self-consistent scheme. Our procedure can be summarized with the 
equation 

Eo = (inf (vl>::.|r + te + Vne\^'..)) + E^^^^[p^.J, (37) 
where E^^^[p,^t^ ] is calculated with the final density resulting from the minimization in the 



first term on the right-hand-side of Eq. (37). 



To carry out the minimization with respect to the potential v^{r) in Eq. (37) we proceeded 
as follows. We parametrized the potential f^(r) with a simple two-parameter form, by 
adding to the physical external potential a gaussian centered on each atom, ce""^^^. The 
minimization of the expectation {'$'^^\T + Vee + Vne\'^vt^) with respect to the two parameters c 
and 7 is done by calculating at each step full- CI wavef unctions ^(^^ for the hamiltonian with 
electron-electron interaction erf { fir 12) /r 12 and external potential f^(r). All calculations were 
done at the cc-V5Z basis-set level. We also produced with MOLPRO ^50j full CI reference 
results for the physical hamiltonian, for comparison. Our simple parametrization of the 
potential f^, containing only two parameters, is enough to yield at /i = the HP energy 
within 0.5 mH, which is the accuracy we sought in this study. This way, we avoid all the 
well-known problems of the OEP method in finite basis set |51j at the price of obtaining 
only an upper bound for our minimization problem (yet, with the reasonable accuracy of 
0.5 mH). 

In Fig. [7] we report the results for E^^^[p] for three different values of the internuclear 
distance R. The dots (•) are the "exact" values of E^^^[p], i.e., the full-CI total energies 
obtained with MOLPRO minus the energies corresponding to the first term on the right- 



hand-side of Eq. (37), inft,M(\I'(^M|T + Vee + V^el^^^p)- The sohd line is -E^n^^M from the 
Overhauser model, and the dashed line is the LDA result, obtained from the parametrization 
of Ref. [2H], in which El!_^^[p] for the uniform electron gas has been calculated with Quantum 
Monte Carlo methods. We see from this figure that when the system is still dominated by 
dynamical correlation, as in the R = 1.4 and the R = 2 cases, the Overhauser model yields, 
even at p = (i.e. for pure KS DFT), correlation energies with errors of ~ 5 mH (while 
LDA is off by ~ 60 mH), which reduce to 1 mH at p = 0.5 (where the LDA error is still 
~ 10 mH). We focuse here on the value p ~ 0.5 since it is the one commonly used in practical 
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FIG. 7: The short-range correlation energy for range-separated multideterminant DFT as a function 
of the cutoff parameter /x for the H2 molecule at three different values of the internuclear distance 



R. Dots (•) are "exact" values (see text in Sec. VI), solid lines are the results from the Overhuaser 



model, and the dashed lines are the LDA values. 
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applications [9l [101 El [I2l [ISl [HI [15]. When the system starts to be dominated by static 
correlation, as in the R = 4 case, the Overhauser model with the simple screened potential of 



Eq. (21) gives, at /i = 0, errors very close to those of LDA (~ 20 mH), which are still of the 
order of ~ 10 mH a.t ^ = 0.5. As the molecule approaches the dissociation limit, R oo, the 
exact tends to the limiting behavior in which -E^~d[p] = E^^[p], and -EcmdH ~ 

for any /i > 0. This is due to the fact that, as R ^ oo, the long-range only wavefunction 
even at very small fi (i.e., with an infinitesimal interaction), becomes essentially exact 
and equal to the Heitler-London wavefunction, so that the functional should be just equal 
to zero. In this limit, the Overhauser model is wrong for small fi (because, as explained 
it misses the "overscreening" at short range), but, for /i ^ 0, yields E^^^[p] 



in Sec. 



IV 



that go to zero much faster than LDA, as it can be already grasped from the third panel 
of Fig. [7| It is thus still more suitable than LDA to be combined with the range-separated 
multideterminant DFT, but it definitely needs some improvement. 

Notice that the Overhauser model would yield much more accurate results if we were able 



to compute -E'^mdW using in Eq. (20) instead of wf^{ri2) the interaction Wgg(ri2) which 
yields the intracule /^(ri2) associated to the wavefunction This way, we would use the 
information available in to the maximum extent, and we would not have the problems 



associated to the "overscreening" discussed in Sec. IV This possibility will be investigated 
in future work. 



VII. CONCLUSIONS AND PERSPECTIVES 



We have presented a preliminary study of the combination of range-separated multide- 
terminant DFT with the Overhauser model, with an application to the paradigmatic case 
of the H2 molecule. We have first analyzed, by means of very accurate variational wave- 
functions, the failure of the Overhauser model in describing static correlation and we have 
then used it to produce an adapted short-range correlation functional for range-separated 
multideterminant DFT. The results are very good for internuclear distances close to equi- 
librium, and are still encouraging as the molecule is stretched. Indeed, in the dissociation 
limit the exact short-range correlation functional should go to zero for any /i > 0, and the 
Overhauser model yields short-range correlation energies that go to zero faster than LDA 
as fi increases. 
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Future work will address the study of better approximations for the unknown Overhauser 
electron-electron interaction, and the development of a more efficient scheme to combine it 
with range-separated multideterminant DFT. 
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